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Abstract 

We rigorously derive explicit formulae for the pair correlation function of 
the ground state of the free Fermi gas in the thermodynamic limit for general 
geometries of the macroscopic regions occupied by the particles and arbitrary 
dimension. As a consequence we also establish the asymptotic validity of the 
local density approximation for the corresponding exchange energy. At con- 
stant density these formulae are universal and do not depend on the geometry 
of the underlying macroscopic domain. In order to identify the correlation 
effects in the thermodynamic limit, we prove a local Weyl law for the spec- 
tral asymptotics of the Laplacian for certain quantum observables which are 
themselves dependent on a small parameter under very general boundary con- 
ditions. 
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1 Introduction and main results 

A classical result of Wigner and Seitz provides an explicit formula for the "exchange 
hole" in ground states of the free electron gas in a cubical box for a large number 
of particles in the thermodynamic limit, |21| : Although the quantum mechanical 
Hamiltonian does not contain coupling terms between different particles, due to the 
Pauli principle, i.e., the requirement that the quantum mechanical wave function be 
antisymmetric with respect to particle exchange, the statistical distribution of the 
individual electrons in J7 is not independent. This phenomenon can be measured in 
terms of the autocorrelation function P(x, y) which is given by the difference of the 
two-body density p2ix + f ,a; — |) and a statistically independent superposition of 
the one-body densities pilx + |) and pi{x ~ |), see Section \TJ\ for details. 

If the number N of electrons in f2 = (0, L)^ tends to infinity at constant micro- 
scopic density p = N/L^, the Wigner-Seitz formula gives an asymptotic expression 
for the autocorrelation function: 

p/ p f 3{sin{pF\y\) - pF\y\cos{pF\y\)) 

^"^'^^^ A iPF\y\r 

where pF = (Stt^p)^/^ is the Fermi-momentum of the free electron gas with density 
p (in atomic units). 

Corresponding to the exchange hole in the particle density there is the "exchange 
energy" E^{A) which measures the difference in the Coulomb energy of the true 
ground state in the region A G ^ as compared to a statistically independent particle 




distribution, cf. Section \TJ\ for details. Results of Dirac, Bloch, Slater and Caspar 
(cf. [51 m [ini [H] ) relate this exchange energy to the "local density approximation" 



EAA)^^c^ / pt^'ix)dx (2) 
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with Cx = 

While these results are classical by now and can be found in standard text 
books on density functional theory as, e.g., by Eschrig [B] or Szabo and Ostlund 
[17) . rigorous convergence proofs seem to have been obtained only much later. The 
book of Thirring [18^ and the articles of Bach [1] and Graf and Solovej [12] give 
rigorous estimates on the exchange energy (not only for the free system). However, 
the first complete convergence proof for the pair correlation function under both 
Dirichlet and periodic boundary conditions has been obtained by Friesecke in [5]. 
In this paper the author also gives a thorough analysis of boundary layer effects and 
sharp error estimates. However, it appears that all formal derivations and rigorous 
results depend heavily on the explicit knowledge of the eigenfunctions of the one- 
body Laplacian on fl. In particular, these calculations do not seem to be easily 
extendable to general domains fl. 

This situation is reminiscent of - and as we will see in fact closely related to - the 
problem of finding a general asymptotic law for the distribution of the eigenvalues 
of the Laplacian on a general domain C R". For = (0,L)^ it had been well 
known that the number of eigenvalues Ai < -^2 < • ■ • of the Laplacian with Dirichlet 
or Neumann boundary values asymptotically satisfies 

Xk ^ 62/3^4/3^-2^2/3 

for large fc, in the sense that the ratio between the left and the right hand side 
converges to 1 for A; — > oo. On physical grounds, at the beginning of the last 
century Lorentz and Sommerfeld had conjectured that the same formula holds true 
for general domains O C M'^. So, in particular, the asymptotic distribution of the 
eigenvalues should depend on only through the volume This conjecture was 
proved not much later by Weyl, cf. [inilSO]- In general dimension n, his result reads 

fe^oo k \n\ 

a formula now known as "Weyl's law" . It holds true, e.g., for the Dirichlet Laplacian 
on general domains C as well as for the Neumann Laplacian under suitable 
regularity assumptions, which are in particular satisfied if dfl is Lipschitz. See, 
e.g., the recent article of Netrusov and Safarov [Ml on Weyl's law for very general 
domains. 

Denoting the number of eigenvalues less than or equal to some A > by N{X) and 
noting that the volume of the ball Br of radius r about is \Br\ = ir^^'^r^ /r{^ + 1): 
we see that Weyl's formula can equivalently be written as 

(2^)"7V(A) ~ |r! X B^l (4) 

for A — oo and thus relates the number of eigenvalues < A to the classical phase 
space volume of x B^. 

Now if {uk ) is an orthonormal series of eigenfunctions corresponding to ( A^ ) , then 
obviously A^(A) — J2k<\(''^k,Uk)L^ ■ Correspondingly, a local Weyl law gives, gener- 
ally speaking, the asymptotics of localized quantities of the form J2k<\i'^k, Au^) , 
where A is a suitable operator. Results in this direction have been obtained for 
pseudodifferential operators A of degree for instance when the Uk are eigenfunc- 
tions of Schrodinger operators — A -|- with smooth potentials V under growth 
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assumptions V{x) — )■ oo as |a;| oo, see, e.g., Zelditch's "Szego limit theorems" 
[55] and also compare the lecture notes T of Evans and Zworski. The volume term 
|11 X on the right hand side of (|H) then has to be replaced with the localized 

term 

a{x, ^) dx d^, 

{e+v(x)<V>^} 

where a denotes the symbol of A. 

Starting with the seminal paper of Shnirelman |15] there have been asymptotic 
results even for the individual terms (uk, Auk)L^ for the eigenfunctions Uk of the 
Laplacian on Riemannian manifolds under the assumption that the geodesic flow 
on the unit cotangent bundle of M be ergodic. For Euklidean domains with smooth 
boundaries and under ergodicity assumptions of he correponding billiard system 
such results for the (free) Laplacian with Dirichlet boundary values are obtained 
in the fundamental contribution by Gerard and Leichtnam |10] . In contrast, our 
main focus will lie on the asymptotics of the Cesaro means '^^=i{uk,Auk)L^ 
for general domains f2. Indeed, radially averaged limits of these quantities with 
polyhomogeneous A have been considered previously in [TU], too. Note, however, 
that the operators A to be investigated below do not introduce radial averaging of 
the Uk- Moreover, they will themselves contain a small length scale depending on 
N. Also, our results will hold true even for very rough boundaries d^l. 

More precisely, our first main result will be a local Weyl law which involves 
two small scales. Let ft C M" be a bounded domain. We consider boundary value 
problems for the Laplacian on f2 subject to the boundary condition Bu = such 
that —A has a pure point spectrum consisting of the eigenvalues Ai < A2 < • • • 
with corresponding eigenfunctions Uk that form an orthonormal basis of 

— Aufc — XkUk in fl, Buk = on dfl, (5) 

such that the sequence (Afc) of eigenvalues obeys Weyl's law ([3|). Note that, e.g., 
this condition is satisfied for the Dirichlet Laplacian on domains f2 as well as for 
the Neumann Laplacian on if dfl is sufficiently regular, say Lipschitz (see, e.g., 
[H]). Define 

and the sequence = Xj, so that Weyl's formula (jS)) reads limfc^oo k^^"'hk — 



-1/2 



The quantities hk '■— Xf. — > as fc — > 00 describe the natural scale for 
measuring the oscillatory wave-length in Uk- On the other hand, after rescaling to a 
fixed macroscopic domain 51, the autocorrelation function for a typical ground state 
of the free Fermi gas is given (say, for simplicity, in the spinless case) by 



PN{x,y) = 



N Z^^ "fc + 2 ) V 2 



(7) 



This introduces a second small length-scale N The determination of the sum 
in d?]) is non-trivial due to the subtle interplay of the two small parameters: By 

Weyl's law, |VMfc(x)|| ^"^^"^ | scales like /i^^iV-i/" and thus like (A:/A^)i/", so the 

terms Uk{x ± ^^—^^ — ~) ^^'^ expected to be close to Uk{x) for fc ^ iV and highly 
oscillating for k ^ N, but to yield non- negligible contributions to Pn{x, y) for fc of 
the same order as N, where the summation is truncated. 
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We denote by {v,w) — (v,w)i^2 — J^^vw the L^(R")-inner product on of func- 
tions two V and w. Deferring the precise definition of the pseudodifferential operator 
a"{x, hD) obtained by the semiclassical Weyl quantization of a symbol a = a{x, ^) 
of class S to Section the result reads as follows. 



Theorem 1.1 Suppose 51 C K" is a hounded domain and (uk) is an orthonormal 
sequence of eigenf unctions for the boundary value problem ([5]), extended by outside 
fl, with corresponding eigenvalues Ai < A2 < . . . that satisfy Weyl's law. If a G S, 
then 



N 

lim y (ufc, a'^ix, N-^/''D)uk) = (27r)-" / a{x, dx d£. 



k=l 

for 7 as defined in ([6]) . 

Note that 7" — ^^^g^^ and hence 

(27r)-"|r2 X B^l = 1. (8) 

So as a functional of a, the limiting expression is nothing but the uniform distribu- 
tion on ri X i?^. 

The starting point for our investigation of correlation effects will be a conver- 
gence result for the so-called one-body matrix, which, rescaled to fixed ft, reads 



Q^{x,y)^N-^Y.^Uk[x + ^^^u^[x ^ j . (9) 

Observing that the Wigner transformation relates this expression to averaged quan- 
tum observables, by applying Theorem ll.il we will be able to obtain the following 
convergence result. 

Theorem 1.2 Under the conditions of Theorem ll.ll the one-body matrices Qn 'in 
dn]) converge to 

^ 2"/-r(t + i)j„,,(7|y|) 

\nMy\r/^ ^""^ ' 

strongly in L^(R" x R") and boundedly in measure on compact subsets of il x R", 
where J„/2 denotes the Bessel function of the first kind of order ^ and 7 is defined 
by dH). 

Since 7" = \nxBi \ ' ~ X-Bi(7^^') and the Fourier transform of the characteristic 

function of the unit ball is given by xihiy) = (27r)"/^|?/|^"/^ J„/2(|y|)7 the term on 
the right hand side is nothing but 

Also note that, although our assumptions on the boundary dVl are very weak, we 
do obtain strong L^-convergence on the whole phase space R" x R". As expected, 
away from the boundary the convergence is even stronger. 

As will be detailed in Section I2.1i for a ground state of the free Fermi gas with 
m spin states, which is given as a Slater determinant of the (rescaled) single particle 
wave functions V'l , ■ • ■ , iJ^n € L^{ft x S), the (rescaled) one-body matrix is given by 



Q%{x, y, si, S2) = iV-i y he + si )iP,[x ^2 . (10) 
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Throughout we wih assume that the tpi (extended by outside VL) satisfy the bound- 
ary condition B4>i{-,s) = on dft for the bounded domain fl so that the set of 
eigenvahies (uk) corresponding to ([5]) satisfies the assumptions of Theorem 11.11 
The corresponding convergence resuh reads 

Theorem 1.3 Let si,S2 G S. The resettled one-body matrices QAr(-, •, si, S2) 'in 
(|10P corresponding to determinantal ground states of the N-body free Fermi gas 
converge to 

, 2"/^r(f + i)j„/2(Pf|;/l)^..s. , - 
Q 0E,y,si,s2 = [7^77 — r-fw2 w 

strongly in i^(R" x R") and boundedly in mettsure on compact subsets of ft x M", 
where pf — 7™^^^". 

We note that p_F is the "Fermi momentum" of the original unrescaled problem: By 
Weyl's law, 

pj.=7m-i/"= hm 7V-i/"7[m-iA^Ji/" = lim N-^/'^X]'^, ,. (11) 

N^oQ TV— foo L / J 

As a direct consequence of these theorems, we obtain a generahzation of the 
Wigner-Seitz formula ([T]) to general shapes of domains in arbitrary dimensions, 
which for the sake of clarity we first state for the spinless mere boundary value 
problem. 

Theorem 1.4 Under the conditions of Theorem \l.l[ the autocorrelation functions 
Pn in converge to 

2"-ir^(t + i)j^^,(7|y|) 

n^^v) \mi\v\Y '''^^^^ 

strongly in L^iW^ x R") and boundedly in measure on compact subsets of fl x R". 
Moreover, we obtain strong convergence of the so-called one-body density 

N 

pi,w(x) :=7V-i^|Mfc(a;)|2. (12) 

k=l 

to the uniform distribution on fl. (Vague convergence of the one-body density for 
the Dirichlet eigenfunctions on smooth domains is in fact inherent already in [lOj.) 



Theorem 1.5 Under the conditions of Theorem \1.U the one-body densities pi^jv 
in (jl2p converge to p = strongly in and boundedly in measure on compact 
subsets of fl. 

For general spin systems we obtain the following generalized formula for the 
exchange hole: Let 

(cf. Section 12.11 for a derivation) . 



N 



,si\ipi[x 



-,S2 



(13) 
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Theorem 1.6 The rescaled autocorrelation functions in (|13p corresponding to 
determinantal ground states of the N-body free Fermi gas converge to 

p^, . 2""^r^(t + i)J^/,(Pi.|y|) 

= m\mP.\y\)- ^"^"^ 

strongly in L^(W^ x K") and boundedly in measure on compact subsets of ft x M", 
where pf is the Fermi-momentum from (1111) . 

2 sin r— r cos r 



As J3/2i^) = Y ^ '^3/2'°'^ ; in n = 3 dimensions the limiting expression reduces 
to 

„S( , f Hsm{pF\y\) -PF\y\cosipF\y\)) 

/ 2-\l/3 

on with = I J , p = I f2 1 ^ ^ . For the electron gas with to = 2 we indeed 

recover the Wigner-Seitz formula ([1]). 

The statement for the one-body density in the general Fermi gas is analogous 
to the spinless case. As justified in Section \2A] we let 



N 

pf,^(x):=7V-i^^|^,(x,s)r. (14) 

sG5 i=l 

Theorem 1.7 The one-body-densities pf in (|14[) corresponding to determinantal 
ground states of the N-body free Fermi gas converge to p — strongly in and 
boundedly in measure on compact subsets offl. 

Thanks to the strong convergence of proved in Theorem II. 6[ we can also 
generalize the Dirac-Bloch formula ^ for the exchange energy in the local density 
approximation to arbitrary domains. It turns out that, for fixed m, the constant Cx 
computed for the box (0, L)^ is in fact universal in three dimensions. 

Recall the definitions of P^ and pf from (|T3)) resp. (IT4|) and set 

KNiA) = [ (15) 

Jaxr^ \y\ 



Theorem 1.8 Let n ~ 3. For any A CC 17 

lim El^{A) = hm / (pf^)^^' = c^\n\-^/^\A\, 

where Cx = 3(rT7r—y^^- 

^ ^ 327rm 

For the electron gas with m ~ 2 and Cx — ji^Y^^ we obtain that formula ^ is 
indeed valid for general domains ft. 

Having stated all our main results, we end this section with a brief outline of the 
following chapters. Section [5] contains the basic material on the physics of the free 
Fermi gas and on Weyl quantization, which will be needed throughout this paper. 
In Section|3]we prove the localized Weyl formula of Theorem ll.il for boundary value 
problems. This will be applied in Section |4] to prove Theorem ll.2l which encodes the 
pair correlation effects in the Cesaro averaged system of eigenfunctions. Theorems 
11.41 and ll.Sl are straightforward corollaries. Systems with spins will then be analyzed 
in Section [S] First Theorem 11.31 will be proved by reduction to the spinless case. 
Then similarly as in the spinless case. Theorems 11.61 11.71 and 11.81 will be direct 
consequences of this result. 
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2 The free Fermi gas and Weyl's calculus 



This chapter serves to collect some background material: firstly, on the quantum 
mechanical description of a Fermi gas and, secondly, on Weyl's pseudodifferential 
calculus for quantum observables. 



2.1 Quantum mechanics of the free Fermi gas 

In this paragraph we will briefly recall some basic material on the correlation func- 
tion of the free Fermi gas, in particular, its connection to the one- and two-body 
densities of the associated quantum mechanical wave function. General references 
are, e.g., Hdl]. 

A Fermi gas of N particles in a region il^ C M" is described by an iV-body quan- 
tum mechanical wave function ip : (i7jv x 5)^ — ^ C, (xi, si,X2, S2, ■ ■ ■ , xn , sn) 
ipixi, . . . , Sn) which is L^-normalized, i.e., 

/ \ip{xi,si,X2,S2, ■ ■ ■ ^xn^sn)]"^ dxi ■ ■ ■ dXN 

and satisfies the Pauli principle, i.e., is antisymmetric with respect to particle ex- 
change. Here S — {ai, . . . ,crm} denotes the set of spin variables, e.g., {cri,cr2} — 
{t,4.} for an electron gas with two states of spin: "up" and "down". In addition, 
ijj is subject to suitable boundary conditions, e.g., Dirichlet conditions for particles 
bound to Qjf by an infinitely deep potential well. 

Particular wave functions are given by the so-called Slater determinants 



f4''"Hx,,s,) ■■■ 4''^\xn,sn)\ 



\i^^ff^\x,,s,) ■■■ ^^^"\xn,sn)J 



(16) 



Here the tpf^"'' are N orthonormal single-particle wave functions. If Ai < A2 < . . . 
denotes the eigenvalues of the Laplacian on il^v and (uk) a corresponding set of 
orthonormal eigenfunctions subject to the boundary condition Buk = on 9f2, 
then a particular orthonormal basis of the single-particle wave functions is given by 

Ml (g) (5cti , ■ • • , Wi (Kl Sa„, ,U2® Sai, ■ ■ ■ ,U2® ^(j,„ , U3 (g) Jo-^ , . . . 

with corresponding energies /^i = . . . = /x^ = Ai, Hm+i = ■■ ■ = fJ.2m = A2, 

The wave function ^ is a ground state for the free Fermi gas if it minimizes the 
kinetic energy 



1 ^ f 

-V" / \\'a:,^P{Xl,...,SN)\'^dXl■■■dxN■ 



i=i ses 

It is well known (and not hard to prove) that a Slater determinant is a ground 
state if and only if ip is the Slater determinant of N orthonormal single-particle 
wave functions ipf^'^\ . . . , V'/^"^ corresponding to the N lowest eigenvalues of the 
one-body problem. 

For given ip, the one-body density is obtained by fixing a point x (z flN inte- 
grating over the remaining spatial variables and all of the spin variables: 



Pi 



[x) := N / \'lp{x,Si,X2,S2,...,XN,SN)fdX2---dXN, (17) 
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the two-body density correspondingly by fixing two spatial points x and x and 
integrating over the remaining spatial variables and all of the spin variables: 

:= f^j ^ \lp{x,Si,x',S2,X3,S3, . . . ,xn,sn)\'^ dx3 - ■ ■ dXN- 

(18) 

The pair correlation function p(^«) then measures the difference of the two-body 
density as compared to a statistically independent superposition of the single par- 
ticle densities: 

P(-")(X, ,) := + f . - - f ) - Ip?'' + I) 

If -0 is a Slater determinant as in (116^ it is well known (and easy to show) that 

N 

and 

pf"\x,x') = E^ E i^f"Hx,^i)n^f"^(^',«2)i^ 

se52 l<jj"<Af 



- "^(X, sOV-f (a:, Sl)V'f "Ha;', S2)^f "Hx', S2) 

and, consequently, 



It will be convenient to also introduce the one body matrix expression 
Q^^"\x,y,s,s') -.^ N / ^p(x+^,s,X2,S2,■■.,XN,SN 

S2,....,Sn£S'"'^N 

' ^ ~ |,s',a;2,S2, • ■ ■,xn,sn^ dx2 ■ ■ ■ dx 

so that 



Q("«' (^, ^, s') = E ^^^"""H^' ^) 

i 

Now in order to investigate the convergence properties of Q^^'^^ and hence p(^«) 
for y of order one in the thermodynamic limit — > oo and VLjq — N^/"-Q M", we 
rescale to a fixed macroscopic region il. Noting that the eigenfunctions scale like 
tP^(x) = ^f\x) = Afi/^V'f "HiVi/"x), we arrive at Qf, : E" x R" x 5 x 5 ^ R, 
P| : R" X M" ^ R and pf^ : R" ^ R given by (HU]), dH]) and dH]), respectively 

2.2 Semiclassical Weyl quantization 

For easy reference we collect some background material on the Weyl pseudodifferential 
calculus that will be needed in the sequel. We will restrict to the very basic symbol 
class 

S := 5(R^") := {a E C°°(R2") : ||9"a||Loo(R2„) < oo for all multiindices a} 
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of smooth functions with bounded derivatives, which wiU be sufhcient for our ap- 
phcations. For much more general classes and additional material on semiclassical 
analysis see, e.g., |H [71 [TB]. 

For h > a, (small) positive number, the semiclassical Weyl quantization of 
a symbol a G S is the operator a"{x,hD) : 5(R") acting on Schwartz 

functions defined by 

.-,.,.Z„.„.),.(2„-"£e-<-.H„(^.,.f).(„..., 

a"{x,hD) extends to a bounded operator from L^(IR") into itself and one has 
{a"{x,hD))* — li" [x,hD). In particular, if a is real valued, then a"{x,hD) is 
self- adjoint. 

If a,b Cz S are two symbols, then also their product generates a pseudifferential 
operator and there exists a constant C = C{a, b) (independent of h) such that 

||a"(x, hD) 6*(x, hD) - (abrix, hD)\\L2^L2 < Ch. (19) 

Garding's inequality asserts that there exists a constant C — C{a) such that for all 
sufficiently smaU h> and aU u G L^(IR") 

(u, a"{x, hD)u)L2 > -Ch\\u\\l2, (20) 

whenever a G S" satisfies a > 0, i.e., is real valued and non-negative. 

Now if a G 5(M^") is itself a Schwartz function on phase space, then the associ- 
ated Weyl quantization is given by 

a""{x,hD)u{x) ^ / ka,h{x,y)u{y)dy (21) 
for all u G -L^(R"), where the integral kernel ka.h G 5(]R.^") is given by 

fc,,,(x, y) = (27r)-" £^ e^(--^)-« a (^^, h^^ (22) 

3 Localized spectral asymptotics 

In this section we will prove Theorem 11.11 For the following preparatory lemmas 
we fix a G 5 of the form a{x,(,) = ai(a;)a2(C), where ai G C^i^l) and a2 G C°°(M") 
has bounded derivatives of all orders. 

Lemma 3.1 Let < a < imin{l,7}, m G N. Suppose 771, 772, X G C°°(M") are 
bounded and such that supp?/i C B^^a; supp 772 C M" \ B^^a o,nd suppx CC ft. 
There exist constants (3 (independent of k, N and a) and C (independent of k and 
N ) such that for sufficiently large N : 

(i) a/3N < fc < (1 - al3)N implies 

\\{l^7jjal®mrix,N''/''D)xuk\\ <CNh'^, j = 1,2 

(ii) and fc > (1 + al3)N implies 

\\{l<E)fj,al<E) ViTix, N~^/"D)xuk\\ < CNh^. 
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With a slight abuse of notation we will sometimes simply write r]i for the function 

Proof. If N~^^"£^ £ suppr^i, then by Weyl's law we have 

\hM = \h + o{l)r'k-'/"i\ 

< 1(7 + o(l))-i(l + a/3)-i/"iV-i/"CI 

< (1 + a)-f-\l + a/3)^i/"(7 + a) 

for fc > (1 + a/3)N, N sufficiently large. For /3 > suitably chosen, the right hand 
side of the last term is bounded by 1 — a for any < a < ^ , whence 

\hkC\ < 1 - a. (23) 

Analogously one obtains 

\hkC\ >l + a (24) 

for |7V-i/"^| > 7 - a and al3N < fc < (1 - a/3)7V, N sufhciently large. (The 
necessary size of N may depend on a and /3.) 

If i = 1 and fc > (1 + a/3)N or if i = 2 and a/3N < fc < (1 - al3)N, then 

{Vjarjir{x,N-'^-D)xuk{x) 

= (2^)-"^^^e^(--^)-%(7V-i/"o?7.(A^-'/"e)a (^^,A^-^/"e) x{y)My)dyd^ 

= (2^)-"/ e''-<^,N,m,^AOi\hke-'^^ [ e-'y<gk{x,y)dyd^ 
foreachm e N, where /fc,jv,™,^,,(0 = '''^"^ a'^^^Ill^-?"''^"^^ "^^^''^"^) and 3^(^,2/) = 

ai(^)x(y)"fe(y). 

It follows that {-qjarji)^ {x, N-^/'^ D)xuk{x) = for x <^ 2n - n. For general x, 
\{rjjari,)'"{x,N-^/"D)xuk{x)\ is bounded by 

(27r)-"||A.,Ar,™,,,,|U2||J-^(-/^2A-l)™.gfe(a;,.)|lL2 
= (27r)-"/2||/,,^,„^,,,|U.||(-/i2A- l)™5,(a;,.)||L2, 

where J^j^ denotes Fourier transform with respect to the y- variable only: {J-yv{x, ■))('?) 
4„ e-"' 2'w(x, y) dy. Since by ^ and ([Ml) 



||/fc,JV,m.ij||i2 < sup 



7,,(7V-i/"^),y,(Ar-i/"^) 



<a-^"'\\r]jfj,\\L^N\\a2\\h 



a2(7V-i/"^) 



2 



and moreover Lemma [3.2l below shows that there exists a constant C only depending 
on m, ai and x such that 

||(-/i2A-l)™5,(a:,.)|U2 <C/ir 

for all fc, it thus follows that 

|(r,,ar?0"(x,^-i/"D)XMfe(a;)| < C(a)7V/i™. 

□ 
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Lemma 3.2 Let ip{y) — ^x{y) — '^i(^^)x(y)- There exist constants C — C{m,^) 
such that for all k Cz N and x G R" 

\\{-hi^^ir4>uk\\L^<ch^. 

Proof. By induction on m it follows that 

\a\<m 

for suitable smooth functions g^n^a supported on suppip. But then 

\a\<m, 

with C only depending on ai and x- 

Consider a nested sequence il DD Hq DD Hi DD ... DD suppV" and first 
observe that ||wfe||_L2(-Qp), Vwfe||i2(-Qj') < C. The first inequality is clear. For the 
second estimate consider a cut-off function 6 e C^{ft) with < 9 < 1 and 9 = 1 
on rii. Then 

ll^Vufc||i.(f,. = / 9^\Vuk\^ = - f 9^UkAuk-2 I 9uuV9-\Iuk 
Jn Jn Jn 

< h^^ + 2\\9\/uk\\mn)\\ukV9\\L2in) < + C\\9yuk\\mn) 

and so \\^Uk\\L^ni) ^ W^^'UkH^n) < Ch^^. 

Standard interior regularity estimates (see, e.g., [11]) for Ad^Uk = h^'^d°'uk on 
give 

and so '||9"ufc||2,2(Q|^|) < C* by induction. (Alternatively one could of course 
directly refer to regularity estimates for powers of A.) Inserting this estimate above, 
we find 

|a| <m 

□ 

If < e < ^ is small enough, we may choose smooth and bounded cut-off 
functions < 771, 772, Xi ^ 1 on R" such that 

suppryi C B^^g, supp 772 C M" \ -B^-e and suppxi C {x G : d\st{x,dVL) > e} 

which satisfy 

rji + 772 = 1 on R" and x = 1 on suppoi. 
Viewing x as a multiplication operator we can state the following lemma. 

Lemma 3.3 There exist constants c (independent of e and N ) andC{e) (indepen- 
dent of N) such that 

N 

^(wfe, a''(x, N-'/"D)uk) - tTaceixiviaviVix, N-'/"D)x) 

k=l 

< ceiV + C(£)7Vi-i/". 
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Proof. By (ITO)) . 

||a-(a;,7V-i/"i?)-Xa"(x,iV-i/"i?)xllM^ <C^(£)^"'/" (25) 
for some constant C{e). Noting that 

(where here and in the sequel we sometimes drop the argmnent (a;, N^^/^D) of the 
Weyl quantized operators) and choosing (3 according to Lemma |3 . 1 1 with a — e and 
m = 3n, we obtain 

|(ufc, xa"x«fc) - {uk, ximamTxuk)] < C{e)Nhl^ < C{e)N~' (26) 

for e^A^ < k < {1- ef3)N, N sufficiently large. (Note that < 2j~^''k-^ for 
k > ePN if TV is large enough.) On the other hand, 

\{uu,x{mamrxUk)\<C{e)Nhl- (27) 

for A; > (1 + sl3)N, N sufficiently large, again by Lemma [5711 

Now a", [rjiarjiY' and multiplication by x ^-re bounded operators on with 
operator norm bounded independently of N . From (|25p . ([^5]) and ([27]) we therefore 
obtain 



^{uk,a"uk) - ^(ufc, xiViaViV XUk] 



k=l k=l 



<Cie)N^-'^/" + ceN + C{e)N 

k=N+l 



3n 
k 



with c not depending on e. The claim now follows from h^" < 27 "^"fc ^ for large 
k and thus T,T=n+iK^"' ^ 7"^"iV-2 for large N and the fact that x{viam)'^X 
vanishes on functions supported in M" \ fi. □ 
Since 7710771 e C;?°(R2"), {r]iaT]i)'"{x,N-^/''D) is of trace class and thus also 
An — xiVi'^Vi)'" X is of trace class with 



Anv^ J kN{x,y)v{y)dy Vw € 5(R"), 

where fcjv is the integral kernel 

kN{x,y) = (27r)-"^^e^(--^)-«77?(iV-i/"Ox(^)x(y)« 

(cf. (HH)). As /cat e C^(R^"), the trace of Ajv can be computed by integrating the 
diagonal elements of the kernel. (This is in fact true for more general kernels, cf. 
[3].) Hence, 



trace Ajv = j kN{x,x)dx 

= (27r)-" / 7jl{N-^/''0x^ix)a{x,N'^/''^)d^dx 
= (27r)-"7v/ rjl{Ox'{^}a{x,Od^dx. (28) 
For easy reference we summarize our observations so far in the following lemma. 

Lemma 3.4 For a = ai (E) 02 E S with ai £ C^(J7) one has 

N 

lim N-'^y^{uk,a"{x, N-^/"D)uk) ^ {2Try / a{x,C)dxd^. 



k=l 
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Proof. Lemma 13731 and Equation ([28)) imply that 

N 



lim sup 



(2^)-"/ x^^WiiOc^i^, Odxd( 



< ce. 



By our choice of x and 771 we thus have 

TV 



hm sup 

TV— >-oo 



< ce 



with c independent of e. The claim now follows since e was arbitrary. □ 

In order to prove Theorem II. II it remains to show that the restrictive assumption 
on the form of a in Lemma 13.41 can be dropped. 

Proof of Theorem ll. 1\ Let a G S* and suppose first that a{x,^) is real- valued, non- 
negative and vanishes for x ^ flp = {z € : dist{z,dfl) > p} for some p > 0. 
For given e > approximating a uniformly by symbols of the form a^{x,^) — 
J2i jgz" a{Si, where {tpf) is a partition of unity on M" with supp V^f C 

5i + (0,2(5)", 5 sufficiently small, by Garding's inequality (|20p we find 



e. 



As a consequence, the assertion of Lemma 13.41 remains true for the symbol a. 

Now more generally consider a e S', a > and let p > 0. Choose a cut-off 
function 9 G C^{^) with < 6* < 1 and 6* = 1 on fip. Let a' = (1 - 6* (g) l)a and 
suppose first that a > 0. By Garding's inequality (1^ we obtain 



and therefore 



TV 



\imMN-^y"{uk,{a r{x,N'^/"D)uk)L-^ > 0. 



(29) 



fc=i 



Now setting a — mO (g) 1 e 5 for m = ||a||oo, again by Garding's inequality (^(7)) 
we obtain 

{uk. (m-a' - hYix, N-^'^D)uk)L- > -CN-^'^. 



From Lemma 13.41 and Equation ([8]) we can now infer that 

TV 



hmsupTV-i V(wfe,(a')"(x,iV-i/"I?)ufc)i2 

TV->oo 

TV 

< limsupA^-1 V (to - (Mfc,(a)-(a;,7V'i/"I?)iifc)i. 

/V_4oo ^ 



fc = l 



TV->oo 

= m - (27r)-" / a(a;,e)dxdC 
< (27r)-"m|(rj\r2p) x 
which tends to zero as p — ?> 0. 



(30) 
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On the other hand Lemma [3.41 yields 

N 

im 

N 



N 

\im N~^y"{uk, {a")'" {x,N'^/"D)uk)L2 ^ C^T^y a" {x,^) dx d£_, (31) 



where a" = {9 1) a, which converges to {2tt) "/fjxs o.{x,£_)dxd^ as p — ^ 0. 
Summarizing 123) , and (PT|) we see that indeed a = a' + a" satisfies 



JV 

hm Ar-iy(a^(a;,Ar-i/"D)iife,Ufc)i2 = (27r)-" / a{x,Odxd^. (32) 



For general a g S* we may choose G R such that 3?a + /i, 9a + /i > 0. With 
the help of dH) the above calculations applied to^a + fi and 9a + /Lt give □ 



4 Averaged correlation asymptotics of eigenfunc- 
tions 

We wish to apply the results of the previous section to investigate convergence 
properties of the pair correlation function (IT^ . In order to do so, we will investigate 
the unsquared one-body density matrix expression from 



1n{x, y) -.^ N ^^Uk[x-\ ] Uk[x 

k=l 



Proposition 4.1 Q^v converges in distributions to Q with 



2"/^r(f + i)j„/,(7|y|) 



Proof. With the help of (|2T|) and ([22]) it is not hard to sec that the formula 
defines a tempered distribution Tk^N e 5'(R" x R"). Since for a e 5(R" x R") 



(27r)-" / e-™-^a(z, O^fe ( ^ + ^' " ) ( z - " ^' ) dz dC 



where we have changed coordinates according to z = w = N^^^^{y — x) and 

( = A^-i/"^^ Tfc,Ar is the i°°-function 



Tk.NizX) = (2^)-"-F 



Mfe I 2 H ]Uk\ Z 



(0- 



Denoting by TyQ{x,y) the partial Fourier transform with respect to the second 
variable only, we find 



TV 



Um{x,0 (2^)-"J-,Qjv(x,0 =^"'I]7^fc,iv(x,0■ 
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Now observe that Theorem 11.11 implies that, as iV ^ oo, 

(C/Ar,a) ^ (27r)"" / a{x,^)dxd£, (33) 

and so 

QN{x,y) a{x,y) dx dy ^ {2Tr)" / F~'^UN{x,y) a{x,y) dx dy 
^ {2ttY' J ^ UN{x,y)Fy^ a{x,y)dxdy ^ j ^ XnxB-,{x,y) J^y^a{x,y) dx dy 

Xn{x)T^'^XB.,iy) a{x,y) dxdy ^ {2T:y / xn{x)xB:{y) a{x,y) dx dy, 



which shows that Q n converges in distributions to 

^, , Xn{x)x^^{y) xn{x)x^,{y) 2"/^r(f + l)Jn/2h\y\) , . 

because 7" = and x^M = {27Tr/^\y\-"/^J,,/,{\y\) and xb-, = XbA'/i)- □ 

The second auxihary resuh for the proof of Theorem 11.21 is the fohowing. 

Proposition 4.2 Qm is bounded uniformly on compact subsets of Q x M". 
In order to prove this result it will be convenient to consider the function 



/n(x, y) ■■=Qn [x-\ , y 



(35) 



Lemma 4.3 Let fl' CC ft, R> 0. There exists a constant C = C{R) such that for 
all N sufficiently large 



N 



sup \QN{x,y)\^ <CN~^y^\uk{:i 



k=l 



for all X £ fl' . 

Proof. U X e fl' and N is large enough, then x + 2N^^^^^Bii C fl, so that y i-> 
Qn{x, y) — N^^ Y^k=i ^k{x+N^^^"y)uk(x) is smooth on B2R and we may calculate 

N 

^yQN{x, y) = A™iV-i ^ «fc (x + 7V-i/"2/) u^{x) 



fe=i 



N 



J2 KN-^^^'^^M^ + N-^'^y)m:{x). 



k=l 



Weyl's law guarantees that < Ck^^"^. With z = x + N ^/"y we thus find 



B2R 



KQN{x,y) 



dy<N 



N 



= N 



dz 



J2 Ar^~'"/"wfc(x + N-^/'^y) Wix) 

k=l 
N 

J2xTN-^'-/"Wix) uk{z) 

k=l 

= N-^J2 A2.'»iV-4™/" \uk{x)\^ 
fe=i 

N 

<C(m)7V-i^|ufc(a;)|2. 



AT 



dy 



k=l 
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The claim now follows from interior regularity estimates on Bji C B2R for the 
Laplacian and Sobolev embedding. □ 

Proof of Proposition ^.'^ Let Vl' CC fi, i? > 0. By Lemma 1131 

N 

\Qn{x, 0)|2 < ^ |y,(a;)|2 ^ c\Qn{x, 0)1 

fe=i 

for all X e il' and so \Qn{x, 0)| < C. Now using Lemma again we find 

N 

sup < C7V-i^|ufe(x)|2 = C|Q^(x,0)| < C. 

y'^^'' k=i 

Hence Qn and thus also Qat is bounded uniformly on compact subsets of J7 x M". 

□ 

Proof of Theorem \1.2l By Proposition 14.21 it remains to show that the convergence 
Q in Proposition 14. 1 1 is in fact strong in L^. To this end, by Proposition 14. II 
it is enough to prove that ||(5Ar||L2 — > |1Q||l2. 

Abbreviating un '■= — 2 — ^ ^^'^ changing coordinates according to w — x + y^, 
z = X — UN, so that I det V(„, ^)(x,2/)| = N, we find 

\QN{x,y)f dx dy 

= N^^ / Uk{x + yN)W{x-yN)u:^{x + yN)u,n{x-yN)dxdy 

TV „ N 

= N^'^ ^ / Uk{w)u;;^{w)u^{z)um{z)dwdz ^ N^^ ^ (5^„ = 1. 

fc,m=l'^"^" k,m=l 

The claim now follows since also Q{x,y) = """^^l^^^^''^^ (cf. ^) satisfies 

\Q{x,y)\'^ dxdy = J^J / | £^ (y) | ^ dy = -|^L /" |xs^(2/)|^ = 1 

because ^=7-^ = 1. □ 
Proof of Theorem \1.4\ This is immediate from Theorem 11.21 □ 
Proof of Theorem\r^ With (HH), © and ([351) we can write 

N 

Pi,n{x) = N^^ ^ Uk{x)uj:{x) = (5^(2^, 0) = (3jv(a:^, 0). 
fe=i 

Thanks to Proposition 14.21 it suffices to establish the L^-convergence of pi^N to 
p. By Theorem 11.21 Qm. and thus also Qjy, converges to Q in i^(R^"), where 
Qi^'V) = ^TOT^' and so |g^(.,y)|2dy ^ 4„ \Q{;y)\'dy in But 



inxB^I 

as shown in the proof of Lemma 



}N{x,y) 



dy = Pi,n{x) 



and also 

2 



which concludes the proof. □ 
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5 Correlation effects for the free Fermi gas 

In this final section we will prove our correlation results for the free Fermi gas 
with spin: Theorems 11.31 II. 6[ 11.71 and 11.81 Here additional "open shell effects" 
have to be considered, which in part can be analyzed similarly as in the case of a 
three-dimensional box (0,L)^, cf. |8]. 

Proof of Theorem \1.S\ Let ajy be the largest integer not exceeding ^ such that 
Aq„+i > Aa„. Similarly, define to be the smallest integer which is not smaller 
than ^ and satisfies Aajv+i > By Weyl's law, N = aNiri + bp^ = ajsim — hM with 
< &iv, ^ N . Then ipam+p ^ "a+i ®5p for a = 0, 1, . . . , oat — 1, j3 = 1, 2, ... m 
extends to a unitary mapping of spanj-j/'i : i < ONm} onto itself. For x',x" e Q, 
si, S2 G S it follows that 



'^1pi{x',Si)lpi{x",S2) ^ ^ Ua+i{x')dp{si)u^{x")Sp{s2) + Rn 

i=X a=0 13=1 

'^^Ua{x')u;;{x")Ssis2 + Rn, (36) 



where Rn = Rn{x' , x", si, S2) = I]ila„m+i ^li^', si)'0i(x", S2). 

The L^-norm of the rest term can be estimated as follows. Setting j/at 
and changing variables to w — x + yN, z — x — i/n so that | det V(u,^z)(x, y)\ — N 
we have 



^ ^ y2 \Rn{x + yN,x - yN,si,S2)\'^ dxdy 

N 

_ 'ipiix + yN,si) i>i{x-yN,s2) 

■,]= 

Tpj (x + yN, si) ipj (x - yN,S2) dxdy 

N 



N 



= N^'^ E E / '^t{w,si) %l)j{w,si) %l>i{z,S2) 'ipj{z,S2) dwdz 

N 

and hence 

N~'^\\RNix + yN,x - yN,si,S2)\\L^(M^^) (37) 

for si, S2 G 5. 

Thus, by (IMl), (021) and Theorem [US we obtain 



AT 

lim N~'^'S^ ipi{x + yN,si)tpi{x - yN,S2) 



lim N~^y^ Ua{x + yN)u^{x - yN)Ssiso 

A^— foo — ^ 

Q = l 

2"/^r(f + l)J„/2(m'V"7|y|)'^.... / ^ n 



S1S2 
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strongly in L^iW"^). 

To conclude the proof it remains to show that Qf^{-,-, 81,82) is bounded on 
compact subsets of 17 x M". From the corresponding result for the spinless case 
in Theorem 11.21 and p6p we see that it is sufficient to prove that N^^Rn IS 
bounded. Since span{-0a„_|_i, . . . , iPn} is a subspace of spanjuo, (8) : + 1 < Q; < 
Q-Afjl < /5 < w}, there exists a partial isometry U = {Ui^(^a^p^) of this space onto 
span{V'aiv+ii ■ • ■ 5 V'w} such that xjji — J2{a p) Ui,(a,i3)Ua'^Sp. We therefore also have 
the following pointwise estimate: 

\Rn{x',x'\8i,S2)\ 

i—aN'm-\~l a— ajv + 1 /3 — 1 ct'— ajv + l — 1 

= Y K(^')n E i""(^")n 

\a=ajv + l / \Q'=ajv + l / 

With the help of Theorem 11.51 we now see that 



|iV ^R{x + yN,x - yN,si,82)\ < ^aj^^ Y Wa{x + yN)\'^ + \ua{x - yN)\^ 

a—UN+l 



remains bounded on compact subsets of x M" . □ 

The proofs of Theorems 11.61 11-71 and 11.81 are now immediate: 
Proof of Theorem I J. 61 By Theorem 11.31 

P^r ^ 1 In^^ 2»-ir^(t + l)J,^^,(p^|y|) 

pUx,y)^^-Y\Q'(-^y'^^'^^)\ = m\mPF\y\r — "^"^"^ 



ses^ 



strongly in i^(R^") and boundedly in measure on compact subsets of O x M". □ 

Proof of Theorem \1.7\ The proof follows along the lines of the proof of Theorem 
11.51 It suffices to show that pf ^ — >■ p in L^(R"). To see this, note ffistly that 

Q%{x + yN,y, 81,82) with j/jv = 2 ^ converges to {x,y, 8i, 82) strongly in 
by Theorem 11.31 secondly that 

Y / \QNi^ + yN.y.Si,82)f dy = pf^„{x) 

Sl ,S2 to 

and thirdly that {x,y, 81, 82) = to^^Q(x, m^^/"?/)5siS2 : whence 



Y, j |<3'^(2^'2/''5i7S2)|^ rfy TO ^ j Q{x,m ^^"y,8i,82) 



dy = p- 



□ 

Proof of Theorem \1.8\ Let A CC fi. By Theorem 11.71 Pin converges to p = 
strongly in Ly^^(Vl) for every p e [1, 00) and so 

1 /3 

hm c. / {pl^f\x)dx = :i(^) M-^'MV 

N^oD J ^ \62iTm J 
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On the other hand, converges to boundedly in measure on ^ x _Bi and 
strongly in on A x (M" \ Si) by Theorem ll.6[ which shows that 



^l/n|^|l+l/n n/2\ I I I 

With the help of the Schafheitlin formula (see, e.g., [221 p. 403]), this integral can be 
determined explicitly: t'"^ Jn/2^''') ~ (n^~i)7r 't- ^ 2. We therefore obtain 



TV-i-oo 



1/3 



lim / j — ^ — dx dy 



(n2 - l)7rTOi/"|17|i+i/"l I 



^327rm^ 

for n = 3. □ 
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